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Abstract. We study the nodal set of eigenfunctions of the Laplace operator on the right angled isosceles 
triangle. A local analysis of the nodal pattern provides an algorithm for computing the number Un of 
nodal domains for any eigenfunction. In addition, an exact recursive formula for the number of nodal 
domains is found to reproduce all existing data. Eventually we use the recursion formula to analyse a 
large sequence of nodal counts statistically. Our analysis shows that the distribution of nodal counts 
for this triangular shape has a much richer structure than the known cases of regular separable shapes 
^ or completely irregular shapes. Furthermore we demonstrate that the nodal count sequence contains 

, information about the periodic orbits of the corresponding classical ray dynamics. 
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1. Introduction 



More than 200 years ago Ernst Chladni pioneered the study of standing waves with his experiments on 
^ ■ sound figures for the vibration modes of plates [Ij. The sound figures revealed that one may characterize 

O |. the modes by looking at the nodal set - the lines on the plate which do not take part in the vibration 

and which are visualized in a sound figure. For each mode he drew the nodal pattern and counted the 
number of nodal lines and nodal domains. With his work he did not only lay the foundations of modern 
acoustics but also started a thread in theoretical and mathematical physics which lead to such classic 
results as Sturm's oscillation theorem [2j and which continues to this day. 

The mathematical framework starts with the Laplacian A on a compact Riemannian manifold A4 - for 
the purpose of this paper it will be sufficient to consider dimension two. If the manifold has a boundary 
then Dirichlet boundary conditions will be assumed. One studies the eigenvalue problem 

(1.1) -A(p = X(p (p\dM = o. 

■ The solutions define the discrete spectrum of (non- negative) eigenvalues {Aat}^^! which we assume to 

be ordered < Ai < A2 < The corresponding eigenfunctions will be denoted cpN- A nodal domain 

of the eigenfunction (pN is a. connected region in M. where the sign of (pN does not change. We define 
the nodal count diS the number of nodal domains in (f^. The nodal counts {i^n} form a sequence of 
integer numbers which characterizes the vibration modes cpN on the shape A4. In case of degeneracies 
in the spectrum the nodal count is not uniquely defined. This may be overcome in various ways, e.g. by 
fixing a basis (and an order in each degeneracy class). Some results on nodal counts are valid for any 
choice for the basis of eigenfunctions - a famous example is the classic theorem by Courant [3j which 

■ states UN ^ 

More recently it has been proposed [4J that one may use the nodal count sequence to distinguish 
between i. regular shapes where the Laplacian is separable and the corresponding ray (billiard) dynamics 
is integrable, and ii. irregular shapes where the ray dynamics is completely chaotic (see also 0[6l[30[9j). 
In the regular separable case the nodal set has a checker board pattern with crossing nodal lines. The 
nodal count can easily be found using Sturm's oscillation theorem in both variables. In this case many 
properties of the nodal count sequence can be developed analytically - e.g. the statistical distribution of 
the scaled nodal count = /N can be described by an explicit limiting function P(^). This function 
has some generic universal features: P(^) is an increasing function with support < ^ < ^crit < 1 
where ^crit is a system dependent cut-off. Near the cut-off, for ^ < ^crit the distribution behaves as 

In the irregular case no explicit counting functional is known. In this case the nodal lines generally 
do not have any intersections and counting nodal domains relies on numerical algorithms (such as the 
Hoshen-Kopelman algorithm [lOj) that represent the eigenfunctions on a grid of finite resolution. The 
numerical procedure is reliable if the resolution is high enough to resolve the distance between nodal lines 
near avoided intersections [llj. For high lying eigenvalues Xn the algorithm is time-consuming due to 
the increasing grid-size. The numerical experiments have shown that a limiting distribution P(^) takes 
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the form P(^) = S{(_ — ^) where ^ is a universal constant (i.e. it does not depend on the shape). This 
and other numerical findings have been shown to be consistent with a seminal conjecture by Berry [12j 
which states that the statistics of eigenfunctions for an irregular (chaotic) shape can be modelled by the 
Gaussian random wave (a superposition of planar waves with the same wavelengths, random direction 
and random phase). Bogomolny and Schmit [13j realized that the nodal structure of a two-dimensional 
random wave may be modelled by a parameter- free critical percolation model (see also [MirTS]). With this 
heuristic model they were able to derive an explicit value for ^ (and other features of the nodal set) with 
excellent agreement to all numerical data. One interesting implication of the critical percolation model 
is that nodal lines can be described by SLE which has been checked affirmative in numerical experiments 
[ini HZl [IH] . Meanwhile some features of the nodal count have been proven rigorously for random waves 
on a sphere - these rigorous results imply the (5- type distribution for P(^) (but cannot predict the value 

Another interesting applications of the nodal count that we will touch in this paper are inverse ques- 
tions. Two inverse questions have been discussed to some detail: i. Can one resolve isospectrality by 
looking at the additional information contained in the nodal count [2Q( \2T\ [22] ? ii. Can one count the 
shape of a drum [231 [211 [23 1211? other words, does the sequence of nodal counts (ordered by increasing 
eigenvalues) determine the shape of the manifold Ail we refer to the shape rather than the manifold 
itself as the nodal counts are invariant under scaling of Ai. 

Both inverse questions have been answered affirmative for certain sets of shapes and some cases have been 
proven rigorously [211[2S]. However, most recently the first example of a pair of non-isometric manifolds 
with identical nodal sequences was found [22j. 

In some cases it could be shown that the geometrical information is stored in the nodal sequence in a 
way which is very similar to the way it is stored in spectral functions. For instance, the nodal count 
sequence for regular shapes with a separable Laplacian can be described by a semiclassical trace formula 
[23} [24| [25] . This trace formula is very similar to the known trace formulas for spectral functions - it is 
a sum over periodic orbits (closed ray-trajectories) on the manifold where each term contains geometric 
information about the orbit. It has been shown that this trace formula can be used to count the shape 
of a surface of revolution [25j . 

In the irregular case, the existence of a trace formula is an open question (unpublished numerical exper- 
iments by the authors give some support to the existence of such a formula). 

In the present work we continue the thread of research summarized above and consider the nodal set of 
the eigenfunctions of one particular shape: the right angled isosceles triangle (i.e. the triangle with angles 
45-45-90). While this shape is regular with an integrable ray dynamics, the Laplacian is not separable. 

Our main result is an explicit algorithm for the nodal counts. In contrast to the numerical algorithm 
used for irregular shapes our algorithm is exact and does not rely on a finite resolution representation of 
the wave function. Though the algorithm is specific to this shape, the approach may serve as the first 
step to generalize explicit formulas for nodal counts beyond the separable case where very few results are 
currently available. Furthermore, we conjecture a recursion formula that allows very efficient evaluation 
of nodal counts for high eigenvalues . 

In the remainder of the introduction we will introduce the spectrum and the basis of eigenfunctions for 
the right angled isosceles triangle. In section [2] we will discuss the nodal structure of the eigenfunctions 
and state the nodal count algorithm and the recursion formula as our main results. In section [3] we 
apply the nodal count algorithm to compute the distribution P(^) of scaled nodal counts, and discuss 
the consistency of the observed nodal counts with the existence of a trace formula. 

1.1. Eigenvalues and eigenfunctions of the Laplacian for the right angled isosceles triangle. 

Let V cM? he the right angled isosceles triangle of area 7r^/2. For definiteness we choose the triangle as 

V = {{x,y) e [0,^]2 : y<x} . 

The eigenvalue problem is stated by 

-A(p{x,y) = - {dl^dl) (p{x,y) = X(p{x,y) with (p{x,y)\Q^ = 
The spectrum of eigenvalues is given by 

Am,n = + for m^n ^ N*and m > n 
and the corresponding eigenfunctions 

(1.2) ^m,n(^7 y) = sin(mx) sm{ny) — sm{nx) sm{my) 
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form a complete orthogonal basis. 

We denote the nodal count (the number of nodal domains) for (prn,n{x^y) by Vrn^n- Let us order the 
spectrum in increasing order, written as a sequence {Aat}^^!, such that Aat < Aat+i- Here N = Nm,n is 
an integer function of the integers m and n (we will continue to suppress the reference to m and n) and we 
have used a mild abuse of notation by writing Xn = XNm,n — ^m,n- The spectrum contains degeneracies 
of a number-theoretic flavour. For a ^-fold degenerate eigenvalue Xn = Aat+i = • • • = Aat+^-i we define 
Nm,n by ordering the degenerate values by increasing n. This ordering is arbitrary and has been chosen 
for definiteness - none of our results here would change with a different choice. 

In principle one may also be interested in the nodal patterns of arbitrary eigenfunctions in a degeneracy 
class. Indeed many physical applications may imply that the basis functions (pm,n cannot be regarded as 
typical as soon as one looks at a degenerate eigenvalue. However, in this paper we will focus exclusively 
on the nodal counts of the basis functions (pm,n ~ for two reasons: i. Understanding the nodal patterns 
of arbitrary superpositions of the basis functions (prn,n is a much harder problem which does not follow 
naturally from understanding just the basis; ii. this choice of basis is natural for any computations. 

2. The nodal pattern 

The current section describes the main properties of the nodal pattern of the eigenfunctions (pmn- These 
observations are then used in subsection 12.31 to infer an exact algorithm for counting nodal domains 
in the triangle. Eventually, we propose a very efficient recursion formula for the nodal counts of the 
eigenfunctions (pmn in subsection 12.51 

2.1. A tiling structure of the nodal lines. The eigenvalue problem on the triangle possess some 
symmetry properties which are revealed in the nodal pattern of the eigenfunctions, (frn,n- Specifically, 
there are eigenfunctions whose nodal sets show a tiling structure: 

(1) For m > n with (m + n) mod 2 = 0, the eigenfunction (/?m,n is an antisymmetric function with 
respect to the line y = ir — x. This line is therefore part of the nodal set of (prn^n- The com- 
plementary nodal set decomposes into two isometric patterns, each from either side of the line. 
Each of these two patterns is similar to the nodal set pattern of the eigenfunction (pm',n' with 

= (m + n)/2 and n' = (m — n)/2 (figure [2m a)). 

(2) For m > n with gcd(m, n) = d > 1 the nodal set of the eigenfunction (prn,n consists of cP identical 
nodal patterns. Each of these patterns is contained within a sub triangle and they are tiled 
together to form the complete pattern. Each such sub pattern is similar to the nodal set of the 
eigenfunction (pm',n' for m' = m/d and n' = n/d (figure ETTJb)). 

The observations above follow directly from (jl.2p . 




(a) (b) 
Figure 2.1. Two examples for the tiling cases: {a)(pg^^ and (b)(/?2i,6 



2.2. Characterization of the nodal set. Let us now characterize the nodal set of the eigenfunction 
^m,n' We assume that the nodal set of (prn,n does not have the tiling behaviour described in section [2711 
i.e. gcd(m, n) = (m + n) mod 2 = 1. Otherwise, one may reduce the values of m, n, as described above, 
to a smaller pair m\n\ that does satisfy this condition, and study the nodal set of (pm^^n' within the 
reduced triangle. In particular, it is proved in Lemma [T] in the appendix that for m, n which satisfy the 
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condition above, the nodal lines of the eigenfunction (pm,n do not cross. This observation is used below 
to characterize the nodal set. 

We write the eigenfunction (pm,n as the difference of the following two functions 

^m,ni^^y) = sin(mx) sin(ny), 
^m,n{x^y) = sin(nx) sin(m^). 

Their nodal sets are correspondingly 

These are regular checkerboard patterns whose nodal domains are open rectangles and triangles (figure 
Ela)). 

The intersection N^^^ fl N^^^ is the set of points 

Vm,n = I < j < i < m| U \ 0<j<i<n} 

I m Jin J 

(marked with bold points in figure [2^ . The eigenfunction (pmn vanishes at these points. Hence, nodal 
lines pass through them. In the following we analyse the run of the nodal lines of (fmn between the 
points of the set Vm^n- The union N^^^ U N^^^ divides V into cells shaped as rectangles and triangles 
of various sizes. These cells are the connected components of V\ (A^4,n ^ ^m,n)- The nodal set of (pmn 
is contained within the cells in which (p^^ ^ and (p'^ ^ have the same sign. These cells are interlacing in 
the checkerboard pattern formed by ^ U ^. We call them the shaded cells and they appear so in 
figure [221b). 

The connection between the points in Vm,n by nodal lines can be easily determined by going over the 
shaded cells and distinguishing between the following cases: 

(1) A rectangular cell adjacent to two points of Vm^n- A non self-intersecting nodal line connects 
these two points. This is proved in Lemma [2l An example is shown in figure [23Ja). 

(2) A rectangular cell adjacent to four points of Kn,n- Two nodal lines connect the two pairs of 
vertices in either a horizontal or a vertical non-crossing pattern. One can determine whether 
the pattern is horizontal or vertical by comparing the sign of (prn,n at the middle point of the 
rectangle with the sign of (pm,n at one of the neighbouring cells. This is proved in Lemma O 
This lemma also proves that a non-tiling eigenfunction, iprn,m cannot vanish at the middle point 
of the rectangular cell. An example is shown in figure [23Jb). 

(3) A cell adjacent to a single point of Vm^n- This happen only for a cell which is adjacent to 
the boundary of V. The Vm,n point is then connected to the boundary of P by a simple non 
intersecting nodal line. This is proved in Lemma [3l An example is shown in figure [23Jc). 

(4) A triangular cell which do not contain any point of Kn,n- In this case there is no nodal line which 
passes through this triangle. This is proved in Lemma [3l 
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Figure 2.3. Different cases of connecting Vm,n within the shaded sub-domains (a) A 
rectangle with two points from Vm,n (b) A rectangle with four points from Vm,n (c) A 
rectangle with a single point from Vm,n 





(a) 



(b) 



Figure 2.4. (a) The nodal set pattern of (fg^^- (b) The graph 6^9,4 which is produced 
by the counting algorithm. 



2.3. An algorithm for counting the nodal domains. We now describe an algorithm for counting 
iyrn,m the number of nodal domains of (pm,m based on the observations of the previous section. If the 
values of m, n correspond to an eigenfunction with a tiling behaviour we replace them by their reduced 
values: 

(1) For m > n with gcd(m,n) = (i > 1, set the new values of m,n to be = m/d and n' = n/d. 
Set the number of tiles to be d^. 

(2) For m > n with (m + n) mod 2 = 0, set the new values of m, n to be = (m + n)/2 and 
n' = (m — n)/2. Set the number of tiles to be 2. 

The number of nodal domains i^rn,n for the original values of m, n equals to the number of tiles times 
the number of nodal domains of the reduced values. We now proceed, assuming the values of m, n were 
reduced. We create a graph, Gm,n, whose vertices are Vm,n with an additional anchor vertex, vq, which 
stands for the boundary of the triangle, dV. The edges of the graph would stand for the nodal lines 
which connect the vertices of Vm^n- We go over all shaded cells as described above and for each of them 
add either zero, one or two edges to the graph connecting the relevant vertices. The number of vertices 
in a cell determines their connectivity, as described in the previous sectiorfl- The cells which contain a 
nodal line connected to the boundary dV^ would contribute a single edge to the graph connecting the 
relevant vertex of Vm,n to the vertex vq. Figure [231 demonstrate the graph Gm,n which corresponds to a 
certain nodal set pattern. 

Once the graph Gm,n is constructed, the number of nodal domains, Vm,m is given by the number of 
interior faces of the graph plus one. According to Euler's formula for planar graphs, the number of faces 
of Gm,n equals E {Gm,n) - \Vm,n\^c{Gm,n), whcrc E {Gm,n) is the number of edges of Gm,n and c {Gm,n) 
is the number of its connected components. We therefore get 

l^E{Gm,n) - |Kn,n| +c(G^,n), 

which completes the algorithm once c {Gm,n) is calculated {E {Gm,n) and \ Vm,n\ are known at this stage). 



In addition, sampling of ipm,n might be required in the case of a cell adjacent to four vertices. 
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2.4. Boundary intersections and nodal loops. Above we have discussed the nodal count Vrn,n' We 
now introduce two further quantities which reflect the nodal set structure of the eigenfunction (^rn,n- The 
first is the number of intersections of the nodal set of (prn,n with the boundary, dV, which we denote by 
r]m,n' The second is the number of closed nodal lines which neither touch the boundary nor intersect 
themselves or any other nodal lines. We call those nodal loops, and denote their number by Im,n- In the 
case where (pm,n does not have a tiling structure, each nodal line is either a loop or a segment connected 
to the boundary at two points. Hence, the connection between the quantities defined above (in the 
non-tiling case) is given by the following formula 

As an example, in figure 12.41 one can count 779^4 = 10 and Ig^4 = 4. The algorithm described in the 
preceding section can be used to count r]m,n and Im,n- 

(1) The number of nodal loops, Im,m is given as the number of connected components of the graph 
Gm,n minus one. 

(2) The number of nodal intersections, r]m,n, equals twice the number of independent cycles of the 
G component which contains vq. 

It was shown recently ([28j) that the number of boundary intersections of the nodal set of ^mn in the 
non-tiling case is given by 

(2.2) r]rn,n = m^n-3. 

Combining this with (j2.ip indicates that any formula for the nodal loop count Im,n would yield a 
formula for the nodal count i'm,n and vice versa. 



2.5. A recursive formula for the nodal loop count. In subsection 12. 3[ we have described an exact 
algorithm that allowed us to compute the nodal loops count. By direct inspection of tables of evaluated 
loop counts we have noticed strong correlations between the counts of different eigenfunctions. An 
extensive analysis of such tables allowed us to infer a recursive formula that we will now describe. Apart 
from regenerating all data that we looked at explicitly, we have checked that the empirical formula 
correctly predicts all loop counts for the first 100,000 non-tiling eigenfunctions (this assures agreement 



of the nodal counts at least up to = 246062, i.e. for all 
We propose that the loop count Im,n is given by 



with 



< 628325). 



Im,n = I (^n, i (m - n - 1) , 

where the 3 parameter function lin^k^l) is defined by the following recursive formula 
(2.3) 

n = lor/c = 

{Ik + (2/ + 1) k'^) + / (nmod {2k + 1) , A:, /) 2A: + 1< n 

21 + 1) (n^ -n) ^i{n,k mod n, /) 2A: + 1 > 2n 

{2P + - n - 2nA: + A:) + ^ A: + / (2A: - n + 1, n - A: - 1, / + 1) n < 2A: + 1 < 2n. 



/(n. A:, /) := < 



As usual we have assumed that m, n correspond to a non-tiling case (otherwise, the reduction described 
above should be made). 



Remarks. 

(1) Note that the description of ()2.3p in terms of the parameters {n^k) = (n, | (m — n — 1)) is more 
compact than a description in terms of the original parameters m, n. 

(2) If the initial values of parameters, n, k correspond to a non-tiling case, i.e. gcd (n + 2A: + 1, n) = 1, 
then this condition will hold for all recursive applications of the formula. 

(3) One can verify that recursive applications of the formula terminate at some stage. Namely, that 
during the recursive applications we arrive at either n = 1 or A: = 0. 
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3. Applications to the nodal counting sequence 

3.1. The nodal count distribution. Let us now discuss the asymptotic statistics of the number of 
nodal domains in terms of the nodal count distribution. In section 11.11 we have given a definition of the 
nodal count sequence {i^n}n=i- Let be the nodal count of the 7V-th eigenfunction. From Courant's 
nodal domain theorem [3j we know that < N. While the Courant bound is only realized by a finite 
number of eigenfunctions [27 \ one may still expect that the nodal count will grow ^ N with the index 
N. It thus makes sense to introduce the scaled nodal count 

(3-1) {» = ^ 

and ask about the asymptotic behaviour of ^at as A/" oo. The latter has been explored by Blum 
et al. [4J for general two-dimensional billiards in terms of the nodal count distribution in the interval 
A < Aat < A(l + g) for large A. The parameter ^ > defines the width of the interval. The limiting 
distribution is defined as 

(3.2) ^^•^(^)=mL) ^ 

where (5e(x) = e (7r(x^ + e^)) ^ is a regularized delta- function (the limit e ^ will always be implied in 
the sequel) and N{X^g) is the number of eigenfunctions in the interval. The integrated distribution will 
be denoted by 

(3.3) hJO= t Px,gAOde- 

Jo 

As mentioned in the introduction an explicit formula for the limiting distribution 

(3.4) PiO = lim PxAO 

A^oo 

can be derived for separable Laplacians using semiclassical methods [4] while for irregular (chaotic) shapes 
Bogomolny's percolation model [13j predicts that the limiting distribution is concentrated at a universal 
value ^ which is consistent with all numerical data available. The right angled isosceles triangle is neither 
an irregular shape (in fact the ray dynamics is integrable) nor are its wave functions separable. The 
proposed recursion formula (j2.3p allows us to find the nodal counts for large sequences of eigenfunctions 
very efficiently on a computer. We calculated the nodal counts for all eigenfunctions with >/A/v < 13000 
(about 66 million eigenfunctions) and extracted the nodal count distributions in various intervals. In the 
remainder of this section we will set ^ = 1 and discuss the numerical results. 

Figure ISTTl reveals that the nodal count distribution Pa,i(0 (with A = 9000^) for the isosceles triangle 
contains a lot of puzzling structure that neither resembles the monotonic behaviour known from separable 
billiards nor the single delta-peak known to describe chaotic billiards. Instead the distribution consists 
of many peaks whose strengths and distances form a visible pattern. Each peak apparently has a further 
substructure. The same structure appears if one only includes wave functions without tiling behaviour 
(or with a specific number of tiles). 

Comparing the nodal count distributions Pa,i(0 various values of A gives us some insight into the 
asymptotic behaviour of Pa,i(0- Figure [321 shows how two peaks in the distribution move and change 
shape as A increases: all peaks move to the left and become sharper. The comparison reveals that 
our numerical calculation of P(0 has not converged - in spite of the extensive number of nodal counts 
included we cannot be sure whether a limiting distribution exists. Still it is interesting to note that, in a 
certain sense, the asymptotic behaviour contains some features of chaotic systems. In a chaotic billiard 
one sees a single peak which becomes a delta- function as A ^ oo. For the isosceles triangle we see a large 
number of peaks - and the numerics suggests that each one may converge to a delta-function. Another 
obvious question suggested by the numerics is whether the limiting distribution contains fractal features. 



3.2. The cumulative nodal loop count. We have already observed in section [2l4l that, at least for 
the non-tiling case, the nodal count decomposes to the number of boundary intersections and the nodal 
loop count (j2.ip . The number of boundary intersections for the triangle was already investigated in [28j 
and presented as a trace formula. In this section we thus focus on the nodal loop count. Denoting by in 
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Figure 3.1. Upper panel: the nodal count distribution (histogram) for energies in the 
interval 9000^ < Aat < 2 • 9000^. The colour represent proportion of wave functions 
with no tiling behaviour (light green), with exactly 2 tiles (dark green), with 4 to 9 
tiles (turquoise), with 10 to 99 tiles (blue), with 100 to 999 tiles (violet), with 1000 to 
9999 tiles (grey), and with more than 10000 tiles (red). Lower panel: the corresponding 
integrated nodal count distribution. 
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P(^) 




Figure 3.2. A detail of the nodal count distribution Pa,i(0 ^^lat shows the limiting 
behaviour. The five curves are histograms for A = 1000^ (orange), A = 2000^ (green), 
A = 4000^ (blue), A = 6000^ (red), and A = 9000^ (black). 

the nodal loop count of the n-th eigenfunction, we define two cumulative continuous counting functions: 

Livj 

n=l 

CO 

n=l 

where [A^J denotes the largest integer smaller than kn is the square root of the n-th eigenvalue 
(multiple eigenvalues appear more than once in the sequence {kn}) and 6 {k) is the Heaviside theta 
function. It should be noted that the functions above can be obtained one from the other by use of the 
spectral counting function, N {k) = Xl^i ^ {k — kn), or its inversion. Previous works examined similar 
nodal counting functions for separable drums [23l [24| . It was proved that for simple tori and surfaces 
of revolution the nodal counting function can be presented as a trace formula. The counting function 
was expressed there as a sum of two parts: a smooth (Weyl) term which reflects the global geometrical 
parameters of the drum, and an oscillating term which depends on the lengths of the classical periodic 
orbits on the drum. For example, it was shown in p3l [24] that the smooth part of i^n is 0{N'^), 

and the oscillating term has the form 




where the sum is over the periodic orbits, Lpo is the length of the orbit, apo, ^po are some coefficients, 
which depend on the orbit, and A is the total area of the drum. Results for other separable drums have 
the same form. 
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-^1,2 -^1,3 -^2,3 Li^4 ^2,4 



Figure 3.3. The power spectrum of Cose (k). The lengths of sorae periodic orbits are 
identified on the / axis. 



Having in mind the case of separable drums, we have examined both Q (N)^ and C (k) numerically 
and found that both counting functions have a (numerically) well-defined smooth term and an oscillatory 
term. Like in the case of the separable drums, the smooth term of C (k) was found to be 0{k^) as well. 
Note that the accumulated boundary intersections count Xl^i rj^Q {k — kn) is only 0{k^). Hence, for 
high energy eigenfunctions, most of the nodal domains do not touch the boundary. We have extracted the 
oscillatory part by numerically interpolating the smooth term and then subtracting it from C (k). In order 
to reveal whether periodic orbits contribute in a similar way as in the separable case we evaluated the 
Fourier transform of the oscillatory term Cqsc {k). The result is shown in figure [331 where the transform 
was performed for the interval 

(^62439153, ^62831853) ^ 

(^\/94662 + 83322, ^^100462 + 7688^^ . 

The Fourier transform in figure 13.31 shows clear peaks at positions which correspond to lengths of 
periodic orbits in the triangle. For each value of (p, g) G \ {(0, 0)}, there exists a continuous family of 
orbits of length Lp^q = 27rY^p2ip^ These are orbits that bounce from the bottom cathetus {y = 0) at 
an angle of arctan((7/p). 

The investigation of Q (N) starts similarly by extracting its oscillating part. As can be expected 
from Weyl's formula, the smooth part is 0{N'^). However, the Fourier transform of Qosc (N) should 
be done with respect to a scaled variable rather than TV. For this purpose, we use the Weyl term of 
the counting function, N ^ '^^n^ where A = ^tt^ is the area of V and Xn is the n-th eigenvalue. 
The scaled variable used for the Fourier transform is the square root of the Weyl-estimated eigenvalue, 

q =^ = ^J ^N. Fourier transforming Qosc with respect to reveals a linear combination of delta- 
like peaks. The positions of these peaks reproduce the lengths of some of the periodic orbits mentioned 
above and of some additional ones: 

(1) Isolated orbits that hit the corner (tt, 0) at 45°. The length of such orbits is Ln = V^Trn, where 
n G N is the number of repetitions of the basic orbit. 

(2) Isolated orbits that go along one of the catheti. Their length is Ln = 27m, where n G N is the 
number of repetitions of the basic orbit. 

Figure 13.41 shows the power spectrum of Qosc{q)^ done when analysing Q (N) in the interval N G 
(38877209, 39269906). 

The above numeric investigation suggests a few observations. The clean Fourier transforms of both 
Cose and Qosc indicate on the existence of a trace formula for both. The need to rescale the variable before 
Fourier transforming Qosc, suggests that the source of a trace formula for Q (N) is the trace formula of 
C {k) combined with the inversion of the spectral counting function, N {k). A similar relation between 
the boundary intersections counting functions was revealed in [28j. Another observation is that only the 
continuous families of periodic orbits appear in the Fourier transform of Cose- This is fundamentally 
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Figure 3.4. The power spectrum of Qosc {q)- The lengths of sorae periodic orbits are 
identified on the / axis. 

different from the trace formula of the boundary intersections ([28]) and calls for further investigation. 
We suggest that the isolated periodic orbits which do appear in the Fourier transform of Qosc are caused 
by the spectral inversion. 

4. Summary & Discussion 

This paper investigates the nodal set of the Laplacian eigenfunctions of the right angled isosceles 
triangle. The novelty of the work is the ability to obtain exact results for the nodal count although this 
problem is not separable. The algorithm described in section 12.31 constructs a graph which reflects the 
topology of the nodal set of a given eigenfunction. The graph contains complete and exact information 
about various properties of the nodal set (such as the number of loops and the number of nodal domains) 
which can be calculated straight forwardly. The standard algorithm used for computing the number of 
nodal domains for a known (non-separable) eigenfunction on a drum is the Hoshen-Kopelman algorithm 
|10| . It samples the eigenfunction on a grid of finite resolution. As far as we know all implementations of 
the Hoshen-Kopelman algorithm for nodal counting use a fixed grid and calculate the number of nodal 
domains as an approximation. In principle, one may reduce the error by increasing the resolution of 
the grid near avoided crossing. However, the application of this algorithm assumes a priori that there 
are no nodal intersections. For the special algorithm we provide here we have proven that it gives the 
exact result, even though it samples the eigenfunction more sparsely than the Hoshen-Kopelman. This 
also leads to a somewhat faster running time of our algorithm (for both algorithms the running time is 
proportional to the energy A - however the constant of proportionality is lower for our algorithm). 

Our result may be generalized to other domains where similar algorithms may apply. Our algorithm 
is based on the fact that the eigenfunctions are presented as a linear combination of simple plane waves. 
It is therefore tempting to try and generalize it for other drums with similar property. The equilateral 
triangle is an immediate candidate (see [29j and references within). 

A further, and quite surprising, result is the recursive formula for the number of nodal loops. To our 
knowledge this is the first known exact formula for the nodal count of a non-separable planar manifold 
(for certain eigenfunctions of tori exact formulas have been given in [22j). The formula was found by 
direct inspection of large tables and has been verified for a large bulk of data computationally. An obvious 
challenge is to prove this formula. In particular, the recursive part of the formula resembles the famous 
Euclid algorithm for the greatest common divisor. A further investigation of the mentioned formula 
might therefore expose some new number theoretical properties of the nodal count. 

The recursive formula enables us to compute a large amount of data and to study the statistical prop- 
erties of the nodal count sequence. We have studied this sequence using functions which are commonly 
used in research of nodal domains: the nodal count distribution and the cumulative nodal count. The 
nodal count distribution showed intriguing structure that resembles neither the behaviour known from 
separable billiards nor the one of chaotic billiards. If at all, there is some similarity to the chaotic case. 
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where the hmiting distribution is a single delta function, whereas in our case it contains a large number 
of peaks. 

In our analysis of the cumulative nodal count we found numerical evidence for the existence of a trace 
formula similar to the one recently derived for separable drums \2T, ^24] . An open question is therefore to 
prove the existence of a trace formula in our case, shedding more light on the question 'Can one count 
the shape of a drum?' 
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Appendix A. Proofs of three lemmas 

Lemma 1. Let 

(A.l) ^m,n(^7 y) = sin(mx) ^m{ny) — sin(nx) sin(m?/), 

he an eigenfunction of the Laplacian on V, where m, n obey gcd(m, n) = {m-\-n) mod 2 = 1. Then there 
are no crossings of the nodal set of (pm,n in the interior of V. 

Proof. The necessary conditions for a crossing to happen at a point (x, y) are 

^m,n {X,y) 0, 
V(/?m,n {X,y) 0. 

After some algebraic manipulations the equations above give 
(A.2) 

(A.3) 

(A.4) 

Combining (jA.2p , (|A.3|1 and (jA.4p gives 

cos {nx) cos {ny) 



sin inx) 


sin 


{ny) 


sin {mx) 


sin 


{myY 


tan {nx) 


n 




tan {mx) 


m ' 




tan {ny) 


n 




tan (m^/) 


m' 





cos {mx) cos (m?/) * 
Squaring this and using (|A.2[) allows to conclude that one of the following holds 

sin^ {my) = sin^ {ny) or 
sin^ {my) = sin^ (mx) . 

Assuming sin^ {my) = sin^ {ny) immediately leads to ^ = ±1, which contradicts the assumptions on 
the values of m and n. Assuming sin^ {my) = sin^ {mx) leads to sin^ {ny) = sin^ {nx). We are now 
required to examine several possibilities for the relations of the expressions mx, my, nx, ny. Such an 
examination shows that each possibility will lead to a contradiction with the requirements x^y G (0, tt) 
and the conditions gcd(m, n) = (m + n) mod 2 = 1. □ 

From now on we consider only m, n obeying the non-tiling conditions. Recall the following definitions. 
Let (pmn be an eigenfunction of the form (|A.ip and 



^mn(^) = sin(mx) sin(ny) 
^mn(^) — sin(nx) sin{my). 



Furthermore 



^m,n = {{x,y)eV\xe-NWye -n| , 

A^^ . = ((x,y) gPIxG-NV^G — n) and 
I ' n m J 



m,n 



m.n 



{^(z, j)|0 < j < z < m} U {^(^, J)|0 < j < ^ < n} . 



By JV{(pmn) we denote the nodal set of (pmn- Let Xc C V \ {N^ ^ U A^ ^) be a rectangular shaped cell 
whose boundary is contained in N^^^ U A^ ^ and contains c points from Vm,ni with po being its centre 
point. We also assume that \/{x,y) G Tc : Sign(pl^^{x^y) = Sign (/9^^(x, 



Lemma 2. 



i J^{^mn) consists of a non- self-intersecting line connecting the nodal corners 0/X2. 

ii ±(pmn{po) > and M{(^mn) nZi consists of two separated lines each connecting adjacent nodal 
corners along edges with ^^mn > 0- 
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Figure A.l. Example of a superposition of the nodal pattern of ip^ and (f^ 

Proof. Nodal sets on 2-dimensional manifolds are submanifolds except for a closed set of lower dimension, 
where nodal lines intersect. For an eigenfunction iprnn this singular set is characterized by ^^^^(0) ^ 
(V(^mn)~^(0). The boundary of a rectangle X2 with two points of Vmn intersects the nodal set only at 
those two points. By elementary arguments using the monotonicity of the sin function, the existence of 
nodal lines that do not intersect with the boundary of this rectangle can be ruled out. The nodal set has 
to connect the nodal corners, since nodal lines do not end. We present this argument in detail for one 
specific case and leave the other cases to the reader. We consider the situation of figure [AT] 

Let the rectangle X2 be in this case such that only one symmetry axis of the two nodal domains of 
ifl^^ and (^^^ enters X2. The symmetry axes are the dotted lines and X2 is shaded. The vertices on 
the lower corners belong to the nodal set and the boundary of X2 between those two points belongs to 
a nodal domain of (^rnn with positive sign (assume this for now - for negative sign it would be the same 
argument). Then the upper boundary of X2 belongs to a nodal domain of ^prnn with a negative sign. 
We study now the behaviour of on a vertical line i between the upper and lower boundary - like the 
one displayed in the figure. Note first that the horizontal rectangle in the figure is a nodal domain of 
^'^r^ with positive sign while the vertical rectangle is a nodal domain of (/:?^^ with positive sign. On the 
lower end of ^, (/:)^^ starts equal to zero and grows strictly monotonic on ^ until it reaches the boundary 
at a positive value, ip]^^ starts with a positive value and falls strictly monotonic ending at zero, iprnn 
being the difference of ip}^^ and (p^^ equals zero exactly once on I. This is true for any t. The nodal 
set therefore intersects every i exactly once and therefore has no intersections nor further isolated nodal 
domains. 

In the case of rectangles with 4 points of Vmn there is a line of constant sign of (frnn running through the 
centre, which cannot be intersected by a nodal line. It can be concluded as above that the nodal corners 
are joined by nodal lines within the two remaining components of this rectangle. □ 

Let T C V \ {N^ ^ U ^) be a triangular shaped cell next to the boundary with Sign(/?^^(x, y) = 
Sign y) in T. Let X^ C V \ {N^ ^ U ^) be a rectangular shaped cell next to the boundary with 

^^g^^lnni^^y) = Sign(/?^^(x,^) in XK 
Lemma 3. 

i A triangular cell T contains a nodal line iffT contains a point ofVm.n- J^i^mn) is a nodal 
line connecting this point to the boundary. 

ii IfX^ contains one point of Vm,n then M{^rnn) H is a nodal line connecting this point to the 
boundary. 

iii IfX^ contains two points ofVm,n then M{(pmn) HX^ is a nodal line connecting those two points. 

Proof. In order to understand the run of the nodal set, the nodal pattern is continued beyond the 
hypotenuse by defining the eigenfunction on the whole square to be the continuation of the eigenfunction 
on the triangle. This rectangle can now be treated just as in Lemma [2] with two points of Kn,n on the left 
lower and right upper corner. The resulting nodal line coincides with the hypotenuse. In case there is a 
point of Vm,n on the right lower corner, there is also one in the left upper corner by symmetry and the 
case with four nodal corners from Lemma [2] applies, and shows the existence of a nodal line connecting 
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the right lower corner with the boundary. The other points are proven similarly to the proof of lemma [2] 
by monotonicity of (f^ and (^^. □ 
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